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We consider the Unambiguous State Discrimination (USD) of two mixed quantum states. We 
study the rank and the spectrum of the elements of an optimal USD measurement. This naturally 
leads to a partial fourth reduction theorem. This theorem shows that either the failure probability 
equals its overall lower bound given in terms of the fidelity or a two-dimensional subspace can be split 
off from the original Hilbert space. We then use this partial reduction theorem to derive the optimal 
solution for any two equally probable Geometrically Uniform (GU) states po and pi — UpoU, 
U 2 — 1, in a four-dimensional Hilbert space. This represents a second class of analytical solutions 
^ ■ for USD problems that cannot be reduced to some pure state cases. We apply our result to answer 

two questions that are relevant in implementations of the Bennett and Brassard 1984 quantum key 
distribution protocol using weak coherent states. 
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Quantum State Discrimination is a crucial task in Quantum Information Theory, especially in a communication 
context. Whenever the signal states are non-orthogonal quantum states, perfect discrimination becomes impossible 
and one must resort to various discrimination strategies. One might, for example, consider an error-free discrimination 
of the states. In that case, due to the non-orthogonality of the signal states, the measurement will sometimes fail to 
identify conclusively the signal. The goal therefore is to minimize the probability of inconclusive results, the so-called 
■ failure probability. This strategy is known as Unambiguous State Discrimination (USD). 

The problem of unambiguously discriminating pure states with equal a priori probabilities was solved by Dieks 
[l[, Ivanovic Q, and Peres 0. Later Jaeger and Shimony presented the general solution for two pure states with 
arbitrary a priori probabilities Q. Shortly after this result, Chefl.es proved that only linearly independent pure states 
can be unambiguously discriminated [f|. Chefles and Barnett then provided the optimal failure probability and its 
corresponding optimal measurement for n symmetric states @. Finally Sun et al Q showed that the unambiguous 
discrimination of pure states is a convex optimization problem [e| 0, [13] • This result was later extended to mixed 
states by Eldar [ljj. With respect to USD of mixed states, three reduction theorems related to simple geometrical 
considerations were developed [l2l [l3l [l4j . They allow us to reduce USD problems to simpler ones where a solution 
might be known. Important examples of reducible problems are Unambiguous State Discrimination of two mixed states 
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with one- dimensional kernel [15j, Unambiguous Comparison of two pure states [l 



17lll8l]. Unambiguous Comparison 



of n pure states with equal a priori probabilities and equal and real overlaps UJ|, State Filtering [191 |20 . |2l| and 
Unambiguous Discrimination of two subspaces [22| . The three reduction theorems also define a standard form of USD 
problems. Such a standard form corresponds to the unambiguous discrimination of two density matrices of rank r in 
a 2r-dimensional Hilbert space without trivial orthogonal subspace and without block diagonal structure [HI, [13, [l4| • 
Interestingly, simple conditions to identify 2x2 block diagonal structures were derived in [23[ ■ When a USD problem 
is not of the standard form, it can immediately be reduced to simpler ones. Necessary and sufficient conditions for 
a USD measurement to be optimal were derived by Eldar J24j. In addition lower and upper bounds on the failure 
probability were provided [la, [25| . Recently, a first class of exact solutions was found [l3| . This class corresponds 
to pairs of mixed states such that the lower bound on the failure probability can be reached. A subclass was later 
studied in 26]. 

In this paper we first study the rank of the elements of an optimal USD measurement. We provide a theorem 
that reveals constraints on the rank of the elements associated to a conclusive detection. If we consider two density 
matrices po and p\ together with their respective a priori probabilities rjo & n( i ?7i, these constraints depend on the 



positivity of the two operators po — w \[PqP\\[Pq an d p\ — y \fPip0\fpt- Note that the positivity of these 
two operators was already introduced in [l~3j as a necessary and sufficient condition for the failure probability to reach 
the overall lower bound 2y / 77t^Tr(y / ^/poPi^7po)- 

A corollary to our first theorem can be derived assuming a standard USD problem. If the two operators po — 
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\/VPnPWPo and p x 




i \fPiPo\fpi arc n °t positive semi-definite, then a two-dimensional subspace can be 



split off from the original USD problem. 

The main result of the present paper is a consequence of our first theorems. We give a second class of exact solutions 
for generic USD problems. This class corresponds to any pair of equally probable Geometrically Uniform (GU) states in 
four dimensions. Two GU states are two unitary similar density matrices po and p\ — Upolf^ where the unitary matrix 
U is an involution i.e. U 2 = 1 . We find that only two options exist depending on the positive semi-definiteness of 
the operator po — \J \/PoPi ^/Po ■ For these two cases we provide the optimal failure probability as well as the optimal 
measurement. We then apply our result to answer two relevant questions related to the implementation of the 
Bennett and Brassard 1984 cryptographic protocol. First 'With what probability can an eavesdropper unambiguously 
distinguish the basis of the signal?' and second 'With what probability can an eavesdropper unambiguously determine 
which bit value is sent without being interested in the knowledge of the basis?'. 

This paper is organized as follows. In Section [IT] we derive a theorem about the rank of the elements of an optimal 
USD measurement and give a corollary which takes the form of a reduction theorem. In Section IIIII we present the 
exact solution for unambiguously discriminating two GU states in a four-dimensional Hubert space. In Section UVl we 
consider two examples of practical interest for quantum cryptography. We then conclude in Section fVl 

II. RANK AND SPECTRUM OF THE ELEMENTS OF AN OPTIMAL USD MEASUREMENT 

In Unambiguous State Discrimination the signal states must be identified without error. If those signal states 
are non-orthogonal quantum states, no perfect discrimination is possible and any USD measurement will lead to 
some inconclusive result. The rate of such inconclusive results is called the failure probability. Given a set of 
signal states pi together with their a priori probabilities rji we want to find an optimal measurement that minimizes 
the failure probability. The measurement is a generalized measurement, that is, a set of Hermitian and positive 
semidefinite operators (called a positive operator-valued measure or POVM [27]) {Ek}k that add up to the identity, 
i.e. J^k — 1- Given a set of N signal states, we consider measurements with N + 1 outcomes, where TV outcomes 
identify conclusively one and only one signal state while the last outcome indicates that the identification failed. The 
N + 1 POVM elements are denoted by Ek, k = 1, . . . , iV, and E?, respectively. The probability to obtain the outcome 
E for a state p is given by the trace quantity Tr(Ep). Therefore the conditions for an error- free measurement simply 
is Tv(EkPi) — whenever k ^ i so that only the state pk can trigger the measurement outcome Ek- The failure 
probability Q to be optimized over all possible USD measurements is given by Q — rjiTv^Efpi). We often note Qi 
the partial failure probability r)iTT(E? p^. 

In this paper we consider the unambiguous discrimination of two signal states po and p± with a priori probabilities 
rjo and 771. Consequently our measurement contains three elements {Eq, E\, E-?} which correspond respectively to the 
conclusive detection of po, to the conclusive detection of p\ and to an inconclusive result. For any Hermitian and 
positive semi-definite operator A, we can introduce the notions of support, kernel and square root. The support Sa of 
A is the subspace spanned by the eigenvectors of A (eigenvectors associated with non zero eigenvalues) . Its orthogonal 

complement is called the kernel of A and is denoted Ka- The square root \f~A of A is defined as the unique positive 

1 — 2 

scmi-dcfinitc operator such that V A — A. This square root operator allows us to write decompositions of the form 



for any unitary transformation V. Since the states pi, i = 0, 1 and the POVM elements Ek, k = 0, 1, ? are positive 
scmi-dcfinitc operators, we can introduce their support, kernel, square root, and decompositions of the form given in 
Eqn.CE]). 

We now derive upper bounds on the rank of the elements E and E\ of an optimal USD POVM. The object of 
our first theorem will be to find under which conditions these upper bounds can be reached. To start, we recall that 
the error-free condition Ti[EiPj] = 0, i, j = 0, 1, i 7^ j, implies the orthogonality between the support of Ei and the 
support of pj [I2I ], Therefore SE i C K. Pj and the dimension of the support SE i: or equivalently the rank rg ; of Ei, is 
upper bounded by the dimension of the kernel JC P . : 



In the relevant case of two density matrices without overlapping supports, dim(H) = dim(S po ) + dim(S pi ) and Eqn.([2]) 
and ([3]) simplify to: 



A = A/M f with M = VAV, 



(1) 



r Eo < dim(K. pi ), 
r El < dim(K- Po ). 



(2) 
(3) 



rE < r o, 



(4) 
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tb x < ri, (5) 

where rj denotes the rank of pt. Note for completeness that it has been already shown in [l3[ that the rank of E? is 
upper bounded by the ranks of the two density matrices po and pi ■ 

r E , < min(r ,ri). (6) 

We are now ready to derive our first theorem. For this theorem, we assume two density matrices without overlapping 
supports (this can always be achieved using the reduction mechanism of [Hj]). The theorem states that the POVM 
elements Eq and E\ of an optimal USD measurement have rank ro and r%, respectively, only if the two operators 
Pa — Jl}^ y^/PoPiy/Po and Pi ~ ^/^V ' \/T>iP()\fPi are positive semi-definite. Here comes the precise statement. 

Theorem 1 Constraints on the rank of the two POVM elements Eq and E\ of an optimal USD measurement 
Consider a USD problem defined by two density matrices po and p\ and their respective a priori probabilities rjg and 
rji such that their supports satisfy S pa PI S pi = {0} (Any USD problem of two density matrices can be reduced to such 
a form according to flB]). Consider also an optimal measurement {EQ Pt , E° pt , E° pt } to that problem. Let F$ and 
F\ be the two operators ^ ^/papi^/po and y/y/pipoy/pi- The fidelity F of the two states po and p\ is then given by 
F = Tr(Fo) = Tr(Fi). Let ro and n be the rank of the two density matrices po and p\. 

If the POVM elements E° pt and E° pt have the same rank as the density matrices po and p\, respectively, then 

p a - .f^F > 0, 
Pi Jfji > 0. 

The proof of this theorem relies on Eldar's necessary and sufficient conditions [13] and tools developed in [3]. A 
detailed proof is given in Appendix A. Interestingly Theorem Q] suggests that the two POVM elements Eq and E\ 

have rank ro and r±, respectively, only in a small regime of the ratio around 1. Indeed the positivity of the two 
operators p — y^Fo an d Pi — y^-Fi is only possible when [l3|, [l8| 

F ~ \no ~ Tr(PoPi)' {> 

Note that these boundaries can be made tighter if more knowledge on po an d oi is provided. Such an example of 
tighter bounds is given in Appendix B. Other boundaries were also derived in [251 ]. 

Theorem [T] can also be rephrased as follows: 



Whenever the two operators po — y^\/VPoPiy/Pn an d Pi — y^\^y/PiPoy/Pi are not positive semi-definite, at 

least one of the two POVM elements Ei, i = 0, 1, of an optimal USD measurement docs not have rank r^. In the case 
of a standard form [341, we can then show that: 



If the two operators po — y ^ \^WPopi\/Po and pi — y ^ \J yfpipo \fpi are not positive semi-definite, then a 
two-dimensional subspace can be split off from the original USD problem. 



This corollary of Theorem Q] actually is a fourth but incomplete reduction theorem. 'Reduction theorem' because no 
optimization is required onto that two-dimensional subspace. 'Incomplete', because the existence and the structure 
of this subspace are known but no complete analytical characterization is available yet. The precise statement of this 
Corollary follows. 

Corollary 1 A fourth, incomplete, reduction theorem 

Consider a standard USD problem defined by two density matrices po and p\ and their respective a priori 
probabilities t]q and r\\ (any USD problem of two density matrices can be reduced to such a form according to 
fZl /j. Consider also an optimal measurement {E^ pt , E° pt , E° pt } to that problem. Let Fq and F\ be the two oper- 
ators y^y/poPi^/po and y \fpiPo\fpi- The fidelity F of the two states po and p\ is then given by F = Tr[Fo) = Tr(Fi). 



4 



Po J£F > 



If one of the conditions < * , — is violated, then (9) 

Pi — Jfji > 

there exists a two-dimensional subspace that can be split off from the original Hilbert space. 

This two-dimensional subspace is characterized by a two-dimensional orthonormal basis {|e), |e')} such that either 



E^le) = |e) 



|e) 6 S po and \e') € IC Po and < E° 



e' = e' 



E° pt \e) = E° pt W) = ^i Pt |e) = E° pt \e') = 0, 



or 

f Er\e) = \e) 

\e) E S P1 and \e') E K P1 and I ^ pt |e') = |e') 

( E°r\e) = E°r\e') = E° pt \e) = ET\e') = 0. 

First let us note that this corollary makes the assumption of a standard USD problem. It is in fact not necessary to 
make such a strong assumption to derive the existence of some eigenvector of E? and -Bq/i with eigenvalue 1 because 
Theorem [T] is valid for any pair of density matrices without overlapping supports. Nevertheless Corollary [1] aims to 
be a fourth reduction theorem. It means in particular that, for any given USD problem of two density matrices, we 
would like to apply the four reduction theorems and always end up with the optimal USD measurement. 

The above corollary is a kind of incomplete reduction theorem. A reduction theorem is a theorem that allows us to 
decrease the size of a USD problem by splitting off some subspace onto which no optimization is needed. To have a 
complete reduction theorem here, we would need to characterize |e) and |e') without solving the whole optimization 
problem. So far the existence of |e) and |e') is ensured and their structure is known (they are eigenvectors of E-? 
and -Eo/i with eigenvalue 1). If |e) and |e') were completely characterized in terms of po, Pi, Vo an d T)i, we would 
have a recipe to solve any USD problem. To see that, let us assume that |e) and |e') can be fully characterized and 
start with two generic mixed states. In the following, the exponent (r) denotes the rank of the density matrices after 
reduction. First, we use the first three reduction theorems to bring the problem into its standard form. We then 

check whether the two operators p^ — \ ^jF^- 1 and p^ — \ ^jF[ r ^ are positive semi-definite. If this is the case, 

\ Vo V 7)1 

then we know the optimal failure probability as well as the optimal measurement to perform since this case falls into 
the first class of exact solutions [13]. If at least one of the two operators p^ — J -^tjFq T - > and — i rnTjF^ is not 

positive semi-definite, we can use our last reduction theorem to get rid of a two-dimensional subspace and define the 
new density matrices p£ and p"-[ together with their respective a priori probabilities t]£ 1 ' and r/[ r . At that 

point we check again the positivity of the two operators 1 ' — J 71 ^^ F^ ^ and p{ ^ — y ^jV-i) F[ r ^ of the 
reduced problem. We see here a constructive way to solve any USD problem. If the two operators ' — J ^jjFq ^ 



and p\ r •* — < n ° F[ r never happen to be positive, we end up with only two pure states and can finally find the 
optimal measurement (see Fig. [1]). A detailled proof of Corollary [T] is given in Appendix C. 

So far, there are only two ways to find a complete characterization of the two eigenvectors |e) and |e') involved in 
Corollary [1] The first possibility is to consider a low-dimensional USD problem. The second option is to consider a 
highly symmetric problem. The former case simply is the case of the two pure states where we want to unambiguously 
discriminate two pure states |\l/o) and l^i) with probabilities r/ and rji. Either the operators po — \ jjr \/^/poPi^/pq 



and pi — ^ \/ y/PiPQy/Pi are positive semi-definite or we have |e) £ S Po/1 eigenvector of E-? with eigenvalue 1, and 
|e') 6 IC Po/1 eigenvector of E 1 / with eigenvalue 1. In only two dimensions, there is no freedom and |e) and |e') must 
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FIG. 1: A constructive way to solve any USD problem if |e) and |e'} in Corollary 1 were fully characterized (the exponent '■*"' 
denotes the rank of the density matrices after reduction) 



be l^o/i) an d l^o/i)- The two pure states case is solved extremely elegantly and rapidly thanks to Corollary[TJ If we 
are interested in higher dimensions we must consider a symmetry to give us enough constraints to fully characterize 
|e) and |e'). With the help of the Geometrically Uniform (GU) symmetry it is possible to go up to four dimensions 
and obtain a second class of exact solutions. This solution of generic USD problem is the object of the next section. 



III. SECOND CLASS OF EXACT SOLUTIONS 



Before deriving the main result of this paper, we need to introduce the so-called Geometrically Uniform (GU) states, 
often met in practical applications [35j |. While symmetric states are generated from one generator state and a single 
unitary transformation, GU states are generated from one generator and a group of unitaries [Tl|, HH [H, H^, H3, El[ . 
More precisely, a set of GU states is a set of density matrices {pi}, i = 0, ...,n — 1 such that pi — UipUj where p is 
an arbitrary density matrix called the generator and the set {£^}, i — 0, ...,n — 1 is a set of unitary matrices that 
form an Abelian group. In order not to break the symmetry of the states, we assume that the a priori probabilities 
are all equal to ^. A consequence of the group structure of the set {Ui} is that we can always consider Uq as 
the identity and po as the generator for a given set of GU states. We can therefore always write two GU states 
as po and p\ — Upoll, where U is an involution (i.e. a unitary transformation U such that U 2 = 1 and U' = U) 
with 770 = 771 = 5. Let us finally note that two GU states are two symmetric states since only a single unitary is needed. 

We can now present the main result of this paper, that is, the optimal failure probability for unambiguously 
discriminating two geometrically uniform states in a four dimensional Hubert space and its corresponding optimal 
measurement. 

Theorem 2 Optimal unambiguous discrimination of two geometrically uniform states in four dimensions 
Consider a USD problem defined by two geometrically uniform states po and p\ of rank 2 with equal a priori probabilities 
and spanning a four- dimensional Hilbert space. Let Fq and F\ be the two operators ^/y/popi^/po and */ ^fplpo^fpl- 
The fidelity F of the two states po and p\ is then given by F = Tr(Fo) = Tr(Fi). We denote by and the 
projectors onto the kernel of po and pi. The optimal failure probability Q°P^ for USD then satisfies 



,o P t = {F if po - F > 

1 1 — (2; I po I s ) otherwise, 



(10) 



6 



where P^ U Pj 1 = a|0)(0| - 6|1)(1| and \x) = -^ g ( e iAr s , « l«'l 1 » Vb\0) + y/a\l)). 
The POVM elements that realize these optimal failure probabilities are given in the two respective cases by 

1. E Q = YT X ifp ~{p Q - F ) Vpo^ -1 where £ = po + Pi (11) 
Ex = UE U 
E-> = 1- E - UE U 



2. E 


= \x){x\ 


Ei 


= UE U 


Eh 


= 1-Eq-UE q U 



Proof We consider a USD problem denned by two geometrically uniform states po and p\ = UpaU, U 2 = 1, of 
rank 2, spanning a four-dimensional Hilbert space. Since rank(po) + rank(pi) — rank(po + pi), the two supports do 
not overlap [l4| • Note that this problem is in the standard form as soon as it has no trivial orthogonal directions and 
no block diagonal structures. 

The symmetry between the two states po and p\ allows further simplifications. Actually, Eldar proved in [241 ] that 
the optimal measurement to unambiguously discriminate geometrically uniform states can be chosen geometrically 
uniform, too. Thus the POVM elements have the form: 

Eo , (12) 

E x = UE U, 

E ? =1-E - UE U. 

In addition we know from the first class of exact solutions 13] that, for two density matrices po and p\ with equal a 
priori probabilities and without overlapping supports, 

Q°» = F* {*--%\\- ^ 

This is even a stronger statement than the desired one since we have an equivalence where we only want an implication 
[3(| . Due to the symmetry of the states, the operators po — Fq and pi — F\ share the same spectrum and the above 
conditions reduce to 

gopt = F ^ po _ Fo > (M) 

If po — Fq ^ 0, then Theorem Q] tells us that at least one of the two POVM elements Eo and E\ does not have 
rank 2, the rank of the two density matrices po and p\. Since E\ = UEqU, Eo and E\ have the same rank so 
that if po — Fo ^ then r^ = r^ 1 < 2. If 7\e = = then E-? = 1 and Q = 1, which is clearly not the 
optimal solution as the density matrices are linear independent. Let us now focus on the remaining case r£ — r^ t = 1. 



First of all it is not too difficult to see that an optimal USD measurement such that te q = = 1 and rank(E?) < 2 
must be a projective measurement with rank(E?) = 2 (See Appendix D for details). For such measurements, we have 
Tr(-Eo-Ea) = or simply (x|J7|x) = where -Bo = |a;)(x|. Since |i) lies in K, Pl , this last relation is equivalent to 

(x\P^UPi\x)=0. (15) 

Importantly, the operator P-^UP^ has one positive and one negative eigenvalue whenever po — Fo is not positive 
semi-definite (a detailed proof can be found in Appendix E). In the eigenbasis {|0), |1)} of P^-UP^, we write 

(16) 
and 

P r UPt=(l _°A a,b>0. (17) 
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If we also consider the normalization of \x), we end up with the following system of two equations: 

\a\ 2 a- \p\ 2 b = 
M 2 + l/?| 2 = l ' 

Up to a global phase, it admits a family of solutions parametrized by a relative phase <E> € [0; 2tt[: 



(18) 



{q= . e ** = ,/? = . 1 =}. (19) 
In the eigenbasis {|0), |1)} we can therefore write 

\x) = | V^P* J . (20) 

We then use again the fact that we are interested in the optimal measurement. Note that we already invoked the 
optimality condition when we used Theorem [T] to state that te q = te± < 2 whenever po — Fq ^ 0. So far \x) is a 
valid choice for any USD projective measurement that is GU symmetric. To find the optimal one, we evaluate the 
corresponding success probability Psuccess- Because of the symmetry of the two GU states, Tt(EqPq) = Tr(Eipi) and 
the success probability for unambiguously discriminating the two GU state po and p\ takes the form 

Psuccess = Tt(EoPo) = (x\po\x). (21) 

After a short calculation, we obtain 

P success = — u (&<0|p |0) + o<l|po|l) + 2VriRe((0\p \l)e- 1 *)) . (22) 
a + b V / 

To maximize this success probability, we choose $ such that i?e((0|jOo|l)e~ 1 *) = |(0|po|l)|- $ must therefore be equal 
to Arg((0\p \l)) and finally 

iArg((0|p |l>) 



This completes the proof. 



Let us remark here that as soon as po — Fq f_ the optimal measurement is 2x2 block diagonal even if the two states 
po and Pi are not. In the next section we consider an example of both theoretical and practical interest. In fact, we 
consider the Bennett and Brassard 1984 protocol (BB84 protocol) implemented through weak coherent pulses with a 
strong phase reference. 



IV. APPLICATION OF THE SECOND CLASS OF EXACT SOLUTIONS TO THE BB84 PROTOCOL 

The Bennett and Brassard 1984 cryptographic protocol 132] provides a method to distribute a private key 
between two parties and therefore allow an unconditionally secure communication. We consider in this section the 
implementation of a BB84-type Quantum Key Distribution (QKD) protocol that uses weak coherent pulses with a 
phase reference (33|. In that context, two important questions related to unambiguous state discrimination can be 
addressed. First, 'With what probability can an eavesdropper unambiguously distinguish the basis of the signal?' and 
second 'With what probability can an eavesdropper unambiguously determine which bit value is sent without being 
interested in the knowledge of the basis?' These two questions can be translated in some unambiguous discrimination 
task concerning two geometrically uniform mixed states in a four-dimensional Hilbert space. We answer these 
two questions providing useful insights for further investigations on practical implementations of Quantum Key 
Distribution protocols. Note that the details of all the following calculations can be found in [14} . 
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A. The pairs of mixed states to discriminate 



The implementation of the BB84 considered here makes use of the four quantum optical coherent states 
{|± a ), | ± t a )} [HI], where a coherent state |/3) is characterized by its complex amplitude (3 S C. The first 
question mentioned above refers to the unambiguous discrimination of the two basis mixed states 



Pr = ~(\a)(a\ + \-at)(-a\), 
pi = ^(\ia)(ia\ + \- ia)(-ia\) . 

The second question refers to the unambiguous discrimination of the two bit value mixed states 

1 



Pa 
Pi 



(\a)(a\ + \ia)(ia\) , 
(\-a)(-a\ + \-ia)(-ia\). 



(24) 
(25) 

(26) 
(27) 



Following the idea presented in [331 ] , we can write these four density matrices in a four-dimensional Hilbert space 



Pr = 



C0C2 

c? cic 3 



c\ 
c c 2 c\ 
cic 3 







(28) 



r 2 
1 



-c c 2 




-c c 2 



-C1C3 





-C1C3 



(29) 



and 



Pa = 



r 2 











1 ClC 


r 2 
c i 












r 2 


^0203 


1 C3C 





^c 3 c 2 


r 2 
c 3 



(30) 



Pi 



^CiCo 



-^C2Ci 







-^C0C3 \ 



l-i 



-C2C3 



\ — ^sco 



-^c 3 c 2 



(31) 



where the coefficients q's can be expressed as functions of the mean photon number p, — \a\ 2 of the coherent pulse 



1 

C2 = 7f 



Cl 



C3 



2 y 7 cosh(^) - 


f cos(^), 




! sin(», 


2 -y/ cosh(/i) - 


- cos(m), 


"2 x/sinhl^) 


— sin(/i) 



(32) 
(33) 
(34) 
(35) 



B. USD of the basis mixed states 



The first question refers to the USD of the two density matrices p r and pi. One can actually calculate the spectrum 
of the operator p r — F r and find that 



Spect(p r — F r ) = {max{c , c 2 }, ma^jcj, c 3 }}} 



(36) 
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which is positive for any value of p. Therefore Theorem [2] tells us that Q opt — F for any value of p. One can also 
calculate the fidelity which takes the form F = |cq — c|| + \c 2 — c||. In terms of the mean photon number the optimal 
failure probability (see Fig. |2]) is finally expressed as 

QO P t = e -M (| cos M | + | sin/x|) , Vai. (37) 



Q 
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1 2 3 4 5 ^ 

FIG. 2: Optimal failure probability for USD of the basis mixed states (fj, in photons per pulse) 



C. USD of the bit value mixed states 

The second question refers to the USD of the two density matrices po an d Pi ■ The spectrum of the operator po — Fa 
is now given by 



1 



Spect(p -F )--[l - e _ " ± e^y/l + e 2 ^ - 2e^ cos(2/i) 



(38) 



This spectrum is not always positive (see Fig. [3]). Indeed only in a regime of relatively large p (p > pq w 0.7193 




FIG. 3: Spectrum of the operator pa — Fq for USD of the bit value mixed states (p in photons per pulse) 



photon per pulse), the quantity ^(1 — — e~ 2 ^^J\ + e 2 ^ — 2e^ cos(2p)) is greater than 0. In the regime p > po, 
the positivity of the operator p — F ensures that the optimal failure probability reaches the fidelity bound F, which 
can be expresses as e _/i (See [14| for details). We therefore obtain 



Q opt = e-",V/i> no. (39) 
Note that for p = pq, the POVM elements E$ and E\ have rank 1 since one eigenvalue of po — Fq vanishes. 

In the regime p < po where the operator p — F is not positive semi-definite Theorem [2] tells us that the optimal 
failure probability for unambiguously discriminating the bit value mixed states is 



)opt _ j 



1 

a + b 



(6(0|po|0) + a(l|p |l)+2Va6|(0|po|l>|) 



(40) 
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where P^U P^ = o|0)(0| - 6|1)(1| and \x) = - ? i=( e ^9« IP«l 1 )) % /6|0) + >/a|l». The complete graph for the failure 
probability is shown in Fig. 21 




FIG. 4: Optimal failure probability for USD of the bit value mixed states (p in photons per pulse) 



So far no neat expression for the optimal failure probability is known in terms of p for p < p$. This comes from 
the rather complicated form of the states po and p± and the fact that no analytical expression of Pq and Pi is known. 
In Fig. Owe finally compare the optimal failure probabilities for USD of the basis and the bit value mixed states. 




FIG. 5: Comparison between the optimal failure probabilities for USD of the basis and the bit value mixed states (p in photons 
per pulse) 



V. CONCLUSION 



In this paper we have provided a theorem giving constraints on the rank of the conclusive outcomes of an optimal 
USD measurement. This theorem also implies a partial fourth reduction theorem. This corollary tells us that either 
the failure probability equals the overall lower bound 2y/rjorjiF or a two-dimensional subspace can be split off from the 
original Hilbert space. In fact, this result can be used as a toolbox to go beyond special cases and solve analytically 
USD problems that can not be solved with the help of the three reduction theorems. We have then employed these 
results to derive a second class of exact solutions. This class corresponds to any pair of equally probable geometrically 
uniform states in a four-dimensional Hilbert space. For that class of states we give the optimal failure probability as 
well as the associated optimal measurement. As an application, we have used this result to address two questions 
related to the implementation of the BB84 QKD protocol with weak coherent states. The questions 'With what 
probability can an eavesdropper unambiguously distinguish the basis of the signal?' and 'With what probability 
can an eavesdropper unambiguously determine which bit value is sent without being interested in the knowledge 
of the basis?' are each related to the unambiguous discrimination of a pair of geometrically uniform states in a 
four-dimensional Hilbert space. 
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VI. APPENDICES 



Appendix A: Proof of Theorem [T] 



We now prove Theorem [T] which is concerned with an optimal measurement for unambiguously discriminating two 
mixed states po and p\. We can therefore use the necessary and sufficient conditions derived by Eldar in [24|. We 
recall them here in a language adapted to our calculations. 

A necessary and sufficient condition for a measurement {Ek}, k £ {0,1,?}, to be optimal is the existence of a 
positive semi-definite operator Z such that 



ZE ? = 0, (41) 

E (Z - mPo)E = 0, (42) 

E 1 (Z - rnpi)Bk = 0, (43) 

P^(Z-r] O p o )P^>0, (44) 

P^-(Z-niPi)P^>0, (45) 



where P- is the projector onto the kernel of pi, i = 0, 1. 

The proof of Theorem Q] can be decomposed in five steps. The first step corresponds to the restriction of Eldar's 
necessary and sufficient conditions to the case where the POVM elements Eq and E\ of a USD measurement have 
rank ro and ri, respectively. After some calculations, this will provide us with a first statement 



If the two POVM elements Eq and E\ of an optimal USD measurement have rank ro and r\, respectively, then 

( ZE? = 

3Z>0 such that < P L (Z-?7 1/ 9 1 )P L = . (46) 

{ p^(z- m p )P^ = o 

Therefore, to prove Theorem [TJ we will simplify have to show the equivalence: 

( ZEh=0 { Po - ,/^F >0 

3Z>0 such that { P^{Z - r] lPl )P^ = I V™ . (47) 

{p^(Z- m p a )P^ = [pi-^i^O 

The remaining four steps will show this equivalence. More precisely, the second step makes use of the notion of 
parallel addition. The third step uses the positive semi-definiteness of the unknown operator Z and focuses on the 
equation ZEi = 0. The fourth step is concerned with the optimal failure probability. The fifth and final step uses an 
equivalence already shown in [1 31 ] . 

First step 

In that first step, we only want to prove an implication. We now see the following: 

If the USD POVM is optimal, i.e. Eqn. ([4T|) - (j43)) are fulfilled, and Eq and E\ have rank ro and r\, respectively, 
then the two Hermitian operators P\{Z — rj p )Pj L and Pq~{Z — ?7 1 p 1 )P 1 must vanish. Indeed the situation is the 
following. We consider two positive semi-definite operators A and B, with A having full rank and ABA^ = 0. We can 
see this relation as the form CC^ = with C = AyB. Such an equation CC^ = is equivalent to C = for any ma- 
trix C. Consequently, ABA^ = is equivalent to A\f~B = 0. Finally, since A is full rank, A~ x exists and B must vanish. 

Let us now focus on Eqn. (|4"2"l) where Eq and P^~(Z — riopo)P^- both have support in JC Pl . We assume 
fE a — dim(K,p 1 )(= ro since the supports of po and p\ do not overlap), so that on the subspace /C Pl we can consider 
E as being full rank. We then set A = E and B = P^{Z - r] Q p Q )P^-. Eqn.(|42|) tells us that ABA^ = with A full 
rank on K. Pl thus B = P\(Z — r)opo)Pj L must vanish. 

We can follow the same idea for Eqn. (|4"3")) . We set A = E\ and B = P(f(Z — rjipi)P^. Since we assume 
fEi = dim(K Po )(= r\ since the supports of po and p\ do not overlap), E\ can be considered as being full rank in IC Po 
and therefore P^{Z — r/iPi)P 1 must vanish. 
Therefore we obtain that: 
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If the two POVM elements Eq and E\ of an optimal USD measurement have rank ro and r±, respectively, then 



3Z > such that < 



ZE> = 
P ± (Z- Vl p 1 )P ± =0 

Pt{Z- Vopo)Pi =0 . (48) 
P^(Z- mPl )P^>0 
. P 1 L (Z- Wo )P 1 ^>0 

The above statement finally implies that: 

If the two POVM elements Eq and -Ei of an optimal USD measurement have rank ro and r±, respectively, then 

( ZE? = 

3Z>0 such that < P^{Z - rnpx)Po = • (49) 

It now remains to show the following equivalence to prove Theorem [TJ 

f = f Po - */£>o > 

3Z>0 such that ^ P L (Z-r]i / 9i)P L = { v ^- . (50) 



if - (£ - %p )Pf L = I pi - V f F i > 

Second step 

Since the supports of the two density matrices po and pi do not overlap, we can make use of the notion of parallel 
addition and introduce the full rank operator = (po + pi) -1 Its main property lies in the relation 

PiYT 1 p 3 = piSij , i = 0,l. (51) 

As a consequence we get the interesting equalities 

PoIT 1 = poZ-'P^, (52) 
Ptp^ 1 = Pt- (53) 

Indeed poT,- 1 = p E- 1 (Pi+P^) = po^T W^+Po^r 1 ^ = po^" 1 ^- Moreover, P± = P^l = Pt{ P a + Pi)^ 1 = 
PtpoY,- 1 . 

If we now consider the equation P\(Z — rjopo)Pj L = 0, we can multiply it on the left-hand side by poS -1 and on 
the right-hand by £ _1 po- We then have 

P^r x Pt(z - 7 ?oPo )F 1 ± S-Vo = 0. (54) 

Because of Eqn. (|52p , this implies 

p o Z- 1 (Z- Vo p o )i;- 1 p o =0. (55) 
We can also go in the other direction. Indeed it follows from the previous equation that: 

PtpoZ-\Z - 7 ?0 Po)S-VoA ± = (56) 

which together with Eqn. (fSU)) yields 

Pt{Z-i 1() p Q )Pt = 0. (57) 

Consequently P\(Z — rj p )P^- = and P qY.^ 1 (Z — rjopo)^ -1 Po — are equivalent propositions for two density 
matrices without overlapping supports. The same result is of course true when we swap and 1. Therefore, 
Pq~(Z — 77ipi)P J - = and p\Yr Y (Z — rjipi)^ 1 p\ = are equivalent too. 

We now come back to the equivalence (|50|) . Thanks to the previous discussion, the proposition 

ZE? = 

3Z > such that { P^(Z - r) X px)Pt = (58) 
P 1 ± (Z- Vo p )P 1 ± = 
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can be advantageously replaced by 



or in short 



where we used Eqn. ([51 



ZE 1 = 

3Z>0 such that { piE" 1 ^- 771^1)1;- Vi =0 (59) 
p Q ^ 1 (Z - <q Q po)Yr 1 pQ = 



=Z::U sudl,1,at < p,s- 1 zrv!= w - for ! = o,i ( 6 °) 



Third step 

Since the operator Z is positive, we know there exists an operator Y such that Z = YY^ . We can insert it in the 
relation p^E -1 ZH~ X p; — r)iPi and find, using the decomposition in Eqn.(l), that there exists a unitary transformation 
Wi such that 

W}Y^- 1 p i = ^Wi^p- h i = 0,L (61) 

Moreover E is full rank since po and pi span the total Hilbert space [13]. We can then decompose Z as Z = 
EE-^E^E = poE-^E-Vo + poE^E'Vi + piE-^E'Vo + piTr 1 ZYr 1 p x . This directly yields 

Z = Vopo + V1P1 + V^o^iVpoWq ^1 \ff>i + VVoViVp^WlWo-s/po (62) 

We finally read off Y* as 

Y* = VP^^ + V^VPi (63) 

where Wt = W^Wq. 

We now make use of the relation ZE? = (Eqn. ([3Tjl ) which is equivalent to Y'E? = 0. We can explicitly write 
= with Ft = ^W^^/pH + y/m^/pl and W = W$Wi. Therefore the statement: 

{%Jj] Q 
n-T-^T^n- -r^n- for , — 1 ( 64 ) 
PjZj ZjZj Pi — r\ l p l , ior % — u, 1, 

can be replaced by: 

There exists a unitary transformation W such that 

- vW f ^E, = y/rjly/PiE? . (65) 

Note that this really is an equivalence and it is not difficult to go from (165]) to (l64l) . If a unitary W exists such that 
—^/rjoW^y/poE-? = y/rjiyfpiE? then we can write (y/rjoW^ ^/po + -s/rji-y/pi)E-? = and define the operator Z — YY^ > 
with Y' as y/rja W' y / po + v^T-v/pi • We then immediately obtain that there exists a positive semi-definite operator Z 
such that ZE-> = 0. To recover Eqn. (|64|) it remains to check that pjS _1 .£E _1 pj = rjipi , for i = 0, 1, which can be 
easily done. 

Eldar's conditions together with the assumptions that Ei (i = 0, 1) have rank and po and pi have no overlapping 
supports are now extremely simplified. We shall now find the final equivalence in two more brief steps. 

Fourth step 

We want to prove that: 

There exists a unitary transformation — W such that — y/rjo W' ^fpoE-? = ^/rjl^fpiE-? 
is equivalent to 

Q opt = 2V^TF. (66) 
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It was already proved a very close statement in [TJI : 

Q opt = y/rfiV*VpoEh = y/pJEi (67) 

where V is a unitary transformation coming from a polar decomposition of J~PqJ~P\- 



This known result already implies that: If Q° pt = 2^ffyrf[F, then there exists a unitary transformation W such 
that -y/ffoW^y/ppE-! = yJrji^fpiE-i. 

The other direction is straitforward too. Indeed, from Eldar's equations we notice that Tr(Z) = P°ul C ess since 
Tr(Z) = Tv(ZE ? ) + Tr(ZP ) + Tr(ZPi) = Tt(^E^Z^E^) + Tt{yfE[Z^/E[) = Tr(v^r?oPoV^) + TH\/^7iPi\/#i) 
where we used Eqns. (|4Tj) . (|42]) and (I43|) . Therefore, with the form of Z we found in Eqn. ([62]) (step 3), the optimal 
failure probability is given by 

Q opt = -Vwh(MVP^W^pl)+Tr(^W^)) (68) 
= 2V%mRe[-Tr(Wt^/pI)]. (69) 

The fidelity can be expressed as F = maxy \Tr(W ^fpoy/pi)\ where the maximum is taken over all the unitary 
transformations. This already implies that 

Q opt < 20M7P (70) 

But of course we know that [TH, [Tj| 

Q opt > 2VWUF- (71) 

In other words, if there exists a unitary transformation —W such that — ^fr^W^ ^fp^E-t = ^frji^fpiE-> , then 
Q opt = ^y/WTlF M- 

Fifth step 

The proof is almost done. Indeed we only need that 



Q opt = 2V^TF O { v ^ (72) 

i?i 1 — 

however this equivalence has been already proved in [l3j . This completes the proof. 




Appendix B: Tighter bounds 



It has been already shown in [Tj 
possible when 



flSj ] that the positivity of the two operators p 



'^F Q and Pl 
no u rl 



Tr(P lPo ) 



< 



< 



F 



Tr(PoPi)' 



2a Pi is only 

m 



(73) 



These boundaries were built considering some very general constraints on Qq and Qi [lc 



Vo^(PiPq) <Qq< Vq, 
77iTr(PoPi) < Qi < Vi, 



(74) 
(75) 



where P, denotes the projector onto the support of /o^, i = 0, 1. If more knowledge on the two density matrices po and 
pi is provided, we can obtain stronger constraints on Qo and Qi and therefore tighter boundaries of the regime (|73|) . 
Let us give such an example of stronger constraints on Qo for, say, a POVM having the GU symmetry E\ = UEqII 
where U 2 = 1. Since E C /C Pl , there exists P > in /C Pl such that P/- = Po + P and therefore Pi + P? = Pi + P. 
Moreover the POVM element P? is invariant under U since PP?P = U(l - E - E\)U = (1 - Pi - P ) = P?. Hence, 
P + P ? = P(Pi + E-?)U = P + [/PP. We therefore derive the trace equality 



Tr(P ? ) = 2Tr(P). 



(76) 
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Indeed Tr(E 1 + P ? ) = TV (Pi) + Tr(P) and Tr(P + P?) = Tr(P ) + Tr(P) so that Tr(l) + Tr(P ? ) = Tr(P ) + TV (Pi) + 
2Tr(P). And, for a USD problem in standard form, the equality Tr(l) = Tr(P ) + Tr(Pi) holds. 
We can now consider Qq. Since Pi + P? = Pi + R and Tr(Pipo) = 0, we can write 

Qo = 77oTr(P ?j0o ) (77) 

= T) Tr(E ? po) + J7oTV(Pipo) (78) 
= 77oTr(Pip )+7 ?0 Tr(Pp ). (79) 

The operator P^-p^P^- is positive semi-definite. We can here introduce X m in> its smallest non vanishing eigenvalue. 
It follows that Qo > ?7oTr(Pi ( o ) + ?7oTr(P)A mm . Together with Eqn.j76]) this yields 

Qo > r7oTr(Pi Po ) + ^iTr(P ? ) (80) 

> ?/o Tr(Pi Po ) + ^iTr(P ?(0o ). (81) 
In other words, for any USD POVM such that Pi = UEolI where U is an involution, 

J- ^rain j <-> 

where X m in — Tnin{Spect(P^~ PqPi )} • This represents a tighter lower bound than the one given in Eqn. (|74p . 

Appendix C: Proof of Corollary [Tl 

To prove this corollary we begin with the statement given in Theorem Q] for two density matrices po and pi with 

the same rank r in a 2r-dimensional Hilbert space. If the two operators po — y^Pb and p\ — ^J^F\ are not positive 

semi-definite, Theorem 1 tells us that at least one of the two POVM elements Po and Pi has rank strictly smaller 
than r. Without loss of generality we say that Te < r. Because of the completeness relation P? + Pi + P = 1 
fulfilled by the POVM elements we have on the support S Pa the equality p)P?Po + PoPiPo + PqEqPq = Po- However 
Sex C JC Po so that we are left with 

P P ? P + PoPoPo = Po- (83) 
Furthermore we can consider the spectral decomposition of the Hermitian operator PoPoPo and write 

r-l 

P E P Q =Y J Mh){h\ (84) 
i=l 

r-l 

Po = ^|A l )(A l | + |e)(e| (85) 

i=l 

where |e) completes the r dimensional orthogonal basis of S Po . As a result P?|e) = (1 — Po — Pi)|e) = |e) — — and 
\e) is an eigenvector of P? with eigenvalue 1. Moreover since |e) is eigenvector with eigenvalue 1 the completeness 
relation is already fulfilled onto the subspace spanned by |e). Therefore no optimization is required onto that subspace 
and we can split it off from the original USD problem. If we denote by >S| e } the subspace of S Po spanned by |e), the 
reduced Hilbert space is Tt/S\ e ) and the support S Po looses one dimension. The remaining USD problem to optimize 
concerns p' and p\ originated from the density matrix po and p\. Here p' has rank r — 1 while p[ has rank r. Thanks 
to the second reduction theorem, we can reduce this problem to the one of two density matrices of rank r — 1 in a 
Hilbert space of dimension 2r — 2. Indeed, the subspace ¥i p > Q fl Sp^ is one dimensional and leads to the detection of 
p[ with unit probability fl2| . We call |e') the unit vector spanning this one dimensional subspace. We are left with a 
reduced USD problem in a 2r — 2 dimensional Hilbert space. Importantly, |e') is in JC p ' Q fl Sp^ C K. p ' () = JC Po . Indeed, 
H = S Pa ® JC Po = S p ' o © <S| e ) © K. Po so that, in H' = H/S^ , JC p ' a — JC Po . 

In other words if po — y^Po and p\ — y^j^Pi are not positive semi-definite, then there exists |e) in S Po , eigenvector 

of P? with eigenvalue 1, and |e') in /C Po , eigenvector of Pi with eigenvalue 1. Without the assumption rE a < rg we 

have the general statement that if po — y^Vp) and p\ — A^Fx are not positive then there exists |e) in either S Po or 

S pi , eigenvector of P? with eigenvalue 1 and |e') either in JC po and eigenvector of Pi with eigenvalue 1, or in JC pi and 
eigenvector of Pq with eigenvalue 1. This completes the proof. I 
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Appendix D: Proof of statement about the projective measurement 

We want to show that an optimal USD measurement such that rg = Te 1 = 1 and rank(E?) < 2 is necessarily a 
projective measurement with rank(E'?) = 2. To do so we can introduce the unit vectors \x) G fc pi , \y) € /C po and the 
real numbers x and y in ]0; +co[ (we could in principle restrict x and y to be in ]0; 1] because probabilities are smaller 
than 1) such that 

E = x\x)(x\ > 0, (86) 
Ei = y\y)(y\>0. (87) 

We call S xy the two-dimensional subspace spanned by \x) and \y), P xy the projector onto it and P^r y the projector 
onto its orthogonal complement. From the definition of the subspace S xy and the completeness relation J2k E k = 1 
we have 

P xy E ^ P xy = P xy ( 88 ) 

Therefore rank{P^ y E^P^ y ) = rank(P xy ) — 2 and Ei must be at least of rank 2. However we already know that 
rank(E-?) < 2. Therefore rank{E-?) = 2 and 

Eh = P^. (89) 

We can now consider the subspace S xy only. On that subspace, we have 

E + E 1 = P xy (90) 

that is to say P xy = x\x)(x\ + y\y)(y\. Since P xy is a projector, P X y — P xy ano - ^ follows that 
x 2 |x)(x| + y 2 \y)(y\ + xy(y\x)\y)(x\ + xy(y\x)\y)(x\ — x\x)(x\ + y\y){y\- The off-diagonal terms are equal if 
and only if (y\x) = (a; ^ and y ^ 0) while the diagonal terms are equal if and only if x = y = 1. Therefore our 
POVM is a projective measurement with rank{E-t) = 2. 



Appendix E: Proof of the statement about the spectrum of P 1 UP 1 

Note that the two operators PqUPq and P^UP^ have the same spectrum. Therefore, we need to prove the 
following theorem 

Theorem 3 Spectrum 

If Po — Fo is not positive semi-definite then P^-UP^- has one positive and one negative eigenvalue. 
First, we write the operator p — F as follows: 



Pa - F Q = 




We now introduce the two orthogonal projectors P± = (since U 2 = 1, P± = P± and P± = P±). If PqUPq is 

positive semi-definite so is y/polly/po and we simply have | y/poU y/po\ — y/pEU ^fp^- In that case the operator p — F 
equals • v //0o(l — C/) v / po = 2^/poP-^/po and is positive semi-definite (since it is of the form ABA^ with B > 0). 
Similarely, if PqUPo is negative semi-definite then the operator —^poU^/po is positive semi-definite and we simply 
have | y/poUyfpo | = —^fpoU^/pH. In that case the operator po — Fq equals v / po(l + U)y/po = 2^poP +y /po and is 
positive semi-definite too. The immediate consequence is that if po — Fo is not positive semi-definite, then PqUPq 
has one positive and one negative eigenvalue. 



To complete the proof, we only need the following simple relation between the spectrum of P UP and that of 
P ± UP ± : 
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Lemma 1 For any two GU states, Spec(PolIPo) = — Spec(P^U P^~) ■ 

This last statement can be proved in three steps. The first step corresponds to the derivation of the eigenvalues 
of Po + UPqII from the eigenvalues of PqUPq. The second step corresponds to the derivation of the eigenvalues 
of Po + UPqII from the eigenvalues of P^UPq . In the last step we compare the two different expressions of the 
eigenvalues of Po + UPqII previously obtained. 

First we can see that the eigenvalues and eigenvectors of PqUPq give us the eigenvalues and eigenvectors of Po + 
UPqII. Let us consider \x), an eigenvector of PqUPq with eigenvalue A i.e. PoUPq\x) — X\x). Since Sp up C Sp , we 
also have Po\x) — \x) and therefore PqU\x) — X\x). Let us now consider the operator Po + UPqII together with the 
vectors P±\x). After some calculations we end up with 

(P + UP U)P±\x) = (l±X)P ± \x). (95) 

This means that to the eigenvector \x) with eigenvalue A of the operator PqUPq correspond two eigenvectors P±\x) 
with eigenvalues (1 ± A) of the operator P + UPqII. Following the same idea, one can show that to the eigenvector 
\x ) with eigenvalue A x of the operator P^UP^ correspond two eigenvectors P±\x ± ) with eigenvalues (1 ± A^) of 
the operator P^ + UP^U. 

Next, we link the spectra of the two operators Po + UPqII and P^-JJP^-. One can actually write 

Po + UPoU = (1 - P Q X ) + f7(l - P X )E7 (96) 
= 21- (Pf + UPj-U). 

This allows us to write the spectrum of the operator Po + UPqII not only in terms of Ai and A2, the eigenvalues of 
PqIIPo but also in terms of Af and A^~, the eigenvalues of P^UP^: 

Spec(P + UP U) = {1 ± Ax, 1 ± A 2 } = {1 t A^, 1 T A^}. 
Note here the important swap between the plus and minus signs. 



The final step follows the observation that two eigenvectors of Po + UPqU are in Sp + and two eigenvectors are in 
Sp_ . We know that the two eigenvalues corresponding to the two eigenvectors in Sp + are given not only by 1 + Ai 
and 1 + A2 (Eqn. ([55]) ) but also by 1 — A^ and 1 — A^-. Indeed from Eqns. ([93|) and ([TO]) , we have for an eigenvector 
{x- 1 ) of Pq^PPq 1 with eigenvalue A- 1 : 

(P + PPoP)P+|a; ± > = [21-(P - L + ^P - L tO]P+|ar L ) (97) 
= (l-A^P+l^). (98) 

Since these two pairs of eigenvalues must be identical, we have either 

1 + Ai = 1 - \i (99) 
1 + A 2 = 1 - Xi (100) 



or 



1 + Ai = l-A^- (101) 
1 + X 2 = l-Xi. (102) 

These two cases only differ from their labellings and we finally end up with {Ai, A2} = {— X^, — A^}. In other words, 
we have obtained that 

S P ec{P UP ) = -Spec(Pj-UPj-). 
This completes the proof of the lemma and therefore the proof of the theorem. ■ 
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